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Fermilab Lattice and MILC Collaborations 

We present an update of our calculation of the form factor forZ? — > D*lv at zero recoil, with higher 
statistics and finer lattices. As before, we use the Fermilab action for b and c quarks, the asqtad 
staggered action for light valence quarks, and the MILC ensembles for gluons and light quarks 
(Liischer-Weisz married to 2+1 rooted staggered sea quarks). In this update, we have reduced the 
total uncertainty on J£"(l) from 2.6% to 1.7%. 

At Lattice2010 we presented a still-blinded result, but this writeup includes the unblinded result 
from the September 2010 CKM workshop. 
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1. Introduction 

The Wbc vertex is proportional to the coupling V c b, which is an element of the Cabibbo [p 
Kobayashi-Maskawa [§] (CKM) matrix. Along with the quark masses, it represents the observ- 
able part of the quarks' coupling to the Higgs sector and is, thus, a fundamental part of particle 
physics. The CKM matrix has four free parameters, and it is convenient to choose one of them to 
be (essentially) |V c b|- Consequently, |V c b| appears throughout flavor physics [Q]. 

| V c b | is determined from semileptonic decays B — > XJv, where X c denotes a charmed final 
state. In exclusive decays, X c is a D or D* meson, and the decay amplitudes can be written 

(D(v D )\r»\B(v B )) = VM^[(v B + v D Yh + (w) + {v B -v D yh_(w)}, (1.1) 
(D*(v D ,a)\r^\B(y B )) = VMsJt^e^e^v^hviw), (1.2) 
(D* (v D , a) \^\B(v B )) = iVMB~M^e { v a) {g Vfl (1 + w)h Al (w) - v»a 2 (w) + v^h Ai (w)} } , (1.3) 

where eW is the D* polarization vector, v B and v D (») denote the mesons' 4-velocities, and w = 
v B • v («) is related to the invariant mass of the Iv pair, q 2 = Mg + M^ (t) — 2wM B M D ( t) . The form 
factors h±, hy, and h Ai (i = 1,2,3) enjoy simple heavy-quark limits and are linear combinations of 
the form factors f±, V, and A, used in other semileptonic decays. 
The differential decay distributions are 

dY{B 2 Dm = ^3-l>(^ + ^) 2 (w 2 -l) 3 / 2 |^| 2 |^(w)| 2 , (1.4) 
dT{B ^f m = ^^(M^M^O^w 2 -!) 1 /^^^^^)^^)^ (L5) 

neglecting the charged lepton and neutrino masses. The physical combinations of form factors are 

M B -M D 2^M B M D 2 

H - {W) = l^m 1 ^ 1 (L6) 




H 2 Q (w)+Hl(w)+HHw) 

3^H M ) ' ( } 

where the zero-recoil (w — > 1) limit of J£" is shown. The function %{w) is chosen so that the square 



root in Eq. ( |1 .7| ) collapses to 1 if hy = h Ai = h Al and h Al = 0, as in the heavy-quark limit without 
radiative corrections. Expressions for H±(w), Hq(w), and %{w) can be found in Ref. [J3j] . 

The messy formula for indicates the advantage of the zero-recoil limit for B — > D*lv: 

one must compute only h Al (1), not four functions. In addition, the heavy-quark flavor symmetry is 
larger when vd* = v B , and Luke's theorem applies. For determining |V c b|, the key aspect of Luke's 
theorem is that it helps control systematic errors. In particular, in lattice gauge theories that respect 
heavy-quark symmetry, one can compute h Al (l) with heavy-quark discretization errors that are 
formally A/rag times smaller than those of h Al (w), w ^ 1, or those of 9?(w) even at w = 1. 

Here we focus on B — > D*lv at zero recoil, describing our calculations of J^"(l) = h Al (\). 
Starting in 2001, experimental determinations of |V c b| used a quenched calculation [Q] 

ff/i\_n Ql O+0.024 _|_ a ni £+0.003 +0.000+0.006 n ox 

where the errors stem, respectively, from statistics, matching lattice gauge theory to QCD, lattice- 
spacing dependence, chiral extrapolation, and the quenched approximation. A notable feature of 
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Table 1: Parameters of the MILC ensembles used for heavy-quark physics. Here C denotes the number of 
configurations in each ensemble; (m'^m' ) the asqtad sea-quark masses; m q the asqtad valence masses; K and 
csw the hopping parameter and clover coupling of the heavy quark. Standard nicknames for the lattice 
spacings are noted (a m 0.045 fm is "ultrafine"). Data are being generated on all ensembles for all m q inside 
the {• • •}, but the present analysis uses at most two, namely m q = m'[ and m q = 0Am' s . 



a (fm) Lattice C (am'^am^) 

w0.15 16 3 x48 596 (0.0290,0.0484) 

medium 16 3 x48 640 (0.0194,0.0484) 

coarse 16 3 x 48 631 (0.0097,0.0484) 

20 3 x 48 603 (0.0048,0.0484) 

w0.12 20 3 x64 2052 (0.02,0.05) 

coarse 20 3 x 64 2259 (0.01,0.05) 

20 3 x64 2110 (0.007,0.05) 

24 3 x64 2099 (0.005,0.05) 

w0.09 28 3 x96 1996 (0.0124,0.031) 

fine 28 3 x 96 1946 (0.0062,0.031) 

32 3 x 96 983 (0.00465,0.031) 

40 3 x 96 1015 (0.0031,0.031) 

ps0.06 48 3 xl44 668 (0.0072,0.018) 

superfine 48 3 x 144 668 (0.0036,0.018) 

56 3 xl44 800 (0.0025,0.018) 

64 3 xl44 826 (0.0018,0.018) 

» 0.045 64 3 x 192 860 (0.0028,0.014) 



m q Kb Kc csw 

{0.0484, 0.0453, 
0.0421, 0.0290, 

0.0194,0.0097, 0.0781 0.1218 1.570 
0.0068, 0.0048} 

{0.05, 0.03, 0.0918 0.1259 1.525 

0.0415,0.0349, 0.0901 0.1254 1.531 

0.02,0.01, 0.0901 0.1254 1.530 

0.007,0.005} 0.0901 0.1254 1.530 

{0.031, 0.0261, 0.0982 0.1277 1.473 

0.0124, 0.0979 0.1276 1.476 

0.0093,0.0062, 0.0977 0.1275 1.476 

0.0047,0.0031} 0.0976 0.1275 1.478 

{0.0188, 0.0160, 0.1052 0.1296 1.4276 

0.0072, 0.1052 0.1296 1.4287 

0.0054, 0.0036, 

0.0025,0.0018} 

{0.014, 0.0056, 0.0028} 



Eq. ( 1.8) is that an estimate of the error associated with quenching has been made. Nevertheless, it 
is necessary to incorporate the light- and strange-quark sea. The first calculation with 2+1 flavors 
of sea quarks obtained [||] 

J^(l) = 0.921 ± 0.013 ± 0.008 ± 0.008 ± 0.014 ± 0.003 ± 0.006 ± 0.004, (1.9) 

where, now, the errors stem from statistics, the gD*Dn coupling, chiral extrapolation, discretization 
errors, matching, and two tuning errors. (The catch-phrases for the errors do not have exactly 



the same meaning in Refs. [Q, |5p; for example, the go-DK error in Eq. ( |L8| ) is incorporated into the 
chiral-extrapolation error.) This paper presents an update of the 2+1-flavor calculation, with mostly 
the same ingredients, but with higher statistics and without the second of the tuning errors. 

The new data set is shown in Table [j], based as before on the MILC ensembles |^] with the 
Liischer-Weisz gauge action |0j, with the g 2 N c ^ but not g 2 Nf corrections and the asqtad- 



improved [10] rooted staggered determinant for the sea quarks. For the valence quarks, we use 



the asqtad action for the light quark and the Fermilab interpretation [|11|] of the clover action [ |12| ] 
for the heavy quark. In this report, we use all ensembles in Table [j] with entries for the heavy- 
quark couplings (Kj,, Kc, and csw), except the fine 32 3 x 96 lattice. These data are being generated 
as part of a broad program of heavy-quark physics, including other semileptonic decays JT3| ] and 
neutral-meson mixing and decay constants [|14|]. 
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Improvements to J^"(l) are timely 0], because the values of \Va>\ that follow from inclusive 



decays are in a 2.2a tension with those that follow from Eq. ( |l.9| ) and also from B — > Dlv and (1) 
[|l5|]. The result described below is but one aspect of a resolution of the discrepancy. Others include 
a re-examination of the extrapolation to zero recoil, unquenched lattice-QCD calculations at w ^ 1, 
lattice-QCD calculations by other groups JI^], and the incorporation of higher-order corrections to 
the inclusive decay expressions. 

In Sec. ||, we discuss details of the data and of the data analysis. Because the value of ^"(1) 
has been studied so much in the past, any new analysis could be influenced in subtle human ways. 
To circumvent any such bias, we hide the numerical value of J^"(l) via an offset in the matching 
factor p A cb, explained in Sec. ||. We present our preliminary results, with all sources of uncertainty 
estimated, in Sec. ^ We include the unblinded value here, which was revealed after Lattice 2010 
but before these proceedings. 

2. Data analysis 

As inRef. [|], we aim for the direct double-ratio 

{D*\c Yj Y 5 b\B)(B\bYjY5c\D*) 2 

^ ~ (D*\cnc\D*)(B\b n b\B) ~ ' {2A) 

where the expressions here are all in (continuum) QCD. To this end, we use lattice gauge theory to 
compute the three-point correlation functions 



C B ^ D \t u t s ,t f ) = £(0|^(x ) t / )^ c y^ b (y > f,)^( 0j f f )|o), (2.2) 
x,y 

C B ^ B {t u t s ,t f ) = J^(0\^ B (x,t f )% r ^b(y,ts)^ B (0,ti)\0), (2.3) 
x,y 

C D ^ D \ ti ,t s ,t f ) = £(0|^ D .(x,/ / )* c y 4 * c (y ! /. s )^.(0 ! / i )|0), (2.4) 



x.y 

where and &d* are interpolating operators coupling to the B and D* mesons, and ^cT^^b and 
^cY^iYs^b are improved currents [11,4], [l7|, Qj. Then the lattice ratio 

C B ^ D \0,t,T)C D ^ B {0,tJ) 
Al{) C D *^ D *(0,t,T)C B ^ B (0,t,T) 

should reach a plateau for a range of t, T S> t 3> 1. The relationship between the plateau value of 

1 /2 

Rj[ and h,A x (1) is discussed in Sec. ||. 

With staggered fermions, G B and &jy* couple to both parities, and three-point correlation func- 
tions have four distinct contributions: 

(f^^tj) = £ £(-l) fa (-l)« r -'W-^V^-0 (2.6) 

= A ^ e -Mxt-M Y {T-t) + { _ y) T-t A ^y e -Mxt-M' Y (T-t) + 

^Yf^Y-^t-ttriT-t) + ^f^Y^t-M'AT-t) + . . . (2 ?) 

with time-dependent factors of —1 associated with the states of undesired parity. To reduce the 
magnitude of the oscillating components, we form the combination [Bl] 
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R Al (0,t, T) = \R Ax (0,t,T) + \R Al (0,t,T + 1) + \R Al (0,t + 1, T + 1), (2.8) 

which should tend more quickly to a plateau. The key here is to have tf = T and T + 1. 

The correlation functions and their ratios are analyzed for two light valence quark masses per 
ensemble, namely, m q = m' t and m q = 0.4m' s (or the single m„ when m'j = 0.4m' s ). The choice of a 
fixed m q , here 0.4m' s , for all m' ; matches, by design, our plans for the ultrafine lattice (a ~ 0.045 fm), 
to anchor future analyses even closer to the continuum limit. Typical plateaus are shown in Fig. |T] 
for a coarse, a fine, and a superfine ensemble. As one can see, the plateau in R Al emerges readily, 
and the statistical errors are 1 % or smaller. 



3. Matching, blinding, and discretization effects 

The ratio combination R Al tends to a ratio of matrix elements like M Ax in Eq. ( |2~I| ) but with 
lattice currents. Each current must be multipled by a matching factor Z A or Zy, defined nonpertur- 



batively in Ref. [17]. The lattice ratio R Al must, therefore, be multiplied by a matching ratio 



n 2 

P A cb 



Z^ cb /ZyccZybb. 



(3.1) 



A subset of the collaboration has computed p AC b in the one-loop approximation. The result is very 
close to unity, but the deviation is, or could be, comparable to h Al (1) — 1. Our numerical analysis 
replaces p A <.b with fkindP^cb, where the blinding factor F\jn n d is again close to unity, but known only 
to those engaged in the one-loop calculation. In this way, choices of fitting ranges, etc., cannot be 
influenced by a human desire to (dis)agree with results for J^"(l) already in the literature. 

The HQET-Symanzik formalism used to define the Zj can also be used to control and suppress 
cutoff dependence Jl8|, [l7|]. In the general case, several operators — both corrections to the current 
and insertions of the effective Lagrangian — generate cutoff effects. For details, see, e.g., the dis- 
cussion of Eq. (2.40) in Ref. For zero recoil, vr>* = vb, and the heavy-quark flavor symmetry 



enlarges from U(l) x U(l) to SU(2). The leading discretization errors drop out, and the remainder 



-1/2 

can be found by applying the formulas of Ref. [ 18] to R^ and h Al (1). One finds 



5 l/2 

Pa aK Ai 



h Al (1) + 0(a s aA 2 /m c ) + 0(a,a 2 A 2 ) + 0(q. 
where the last error acknowledges the one-looD calculation of o„ ■ 



(3.2) 



' r~ri r~\i i r- 
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Figure 1: Ratio combination RJ (0,t,T) vs. t with m q = = 0.2m' s . From left to right: the coarse ensemble 
with T = 12 and (am' h am' s ) = (0.01,0.05); the fine ensemble with T= 17 and (am'„am' s ) = (0.0062,0.031); 
the superfine ensemble with T = 24 and (am'^am'^) = (0.0036,0.018). 
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behavior of Fermilab actions provides a reasonable guide to the dependence on rriQa of the correc- 



tions. We see in our data little dependence on the lattice spacing, in accord with Eq. ( p. 4 ) 
4. Preliminary result 



Figure || provides a glimpse into our systematic error analysis, which closely follows Ref. [Q]. 
We use our previous study of heavy-quark-mass dependence to fine-tune a posteriori the hopping 
parameters and to assess the tuning errors. We fit the light-quark mass dependence to one-loop 
chiral perturbation theory, suitably modified for staggered quarks Jl9[] . The cusp is a necessary, 
physical effect that appears because the Dn threshold sinks below the D* mass. 

With the blinding factor in place, we find 



fkmd^(l) = 0.8949 ± 0.0051 ± 0.0088 ± 0.0072 ± 0.0093 ± 0.0030 ± 0.0050, 



(4.1) 



where the errors again stem from statistics, the gD*Dn coupling, chiral extrapolation, discretization 
errors, matching, and tuning and jq,. To show how the errors have been reduced, it helps to scale 
this result to the old central value (Fp is the needed ad hoc factor): 



^(1) = 0.921(13)(8)(8)(14)(3)(6)(4) [g, 
F P &(1) = 0.921 (05) (9) (7) (10) (3) (5) [this work]. 



(4.2) 
(4.3) 



The higher statistics and wider scope of this dataset has reduced the statistical error wither 1 / 2 . The 
quoted heavy-quark discretization error is smaller, because with the superfine data we can move 
beyond pure power counting and combine the (lack of) trend in the data with the detailed theory 



of cutoff effects [|18|]. After Lattice 2010, we continued to examine the heavy-quark discretization 
and ff-tunings errors, reducing them somewhat, and the chiral-extrapolation error, increasing it 
somewhat. For the 2010 Workshop on the CKM Unitarity Triangle, we removed the blinding 



factor, finding [20]: 



[1) = 0.9077(51)(88)(84)(90)(30)(33). 
This result reduces the tension with |V C b| from inclusive decays to 1.6a. 



(4.4) 
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Figure 2: Left: dependence of (1) on the heavy-quark hopping parameters (with data of Ref. 
chiral extrapolation showing only points with m q = mi and a fit to all data. 



)■ Right: 
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